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ABSTRACT 

Acoustic  data  from  the  Arctic  Climate  Observations  using  Underwater  Sound  (ACOUS) 
experiment  are  analyzed  to  determine  the  correlation  between  acoustic  propagation  loss  and  the 
seasonal  variability  of  sea  ice  thickness.  The  objective  of  this  research  is  to  provide  long-term 
synoptic  monitoring  of  sea  ice  thickness,  an  important  global  climate  variable,  using  acoustic 
remote  sensing.  As  part  of  the  ACOUS  program  an  autonomous  acoustic  source  deployed 
northwest  of  Franz  Josef  Land  transmitted  tomographic  signals  at  20.5  Hz  once  every  four  days 
from  October  1998  until  December  1999.  These  signals  were  received  on  a  vertical  array  in  the 
Lincoln  Sea  1250  km  away.  Two  of  the  signals  transmitted  in  April  1999  were  received  on  a 
vertical  array  at  ice  camp  APLIS  in  the  Chukchi  Sea  north  of  Point  Barrow,  Alaska,  at  a  distance 
of  approximately  2720  km  from  the  source.  Temporal  variations  of  the  modal  propagation  loss 
arc  examined.  The  influence  of  ice  parameters,  variations  of  the  sound  speed  profile,  and  mode- 
coupling  effects  on  the  propagation  losses  of  individual  modes  is  studied.  The  experimental 
results  are  compared  to  the  results  of  the  earlier  experiments  and  the  theoretical  prediction  using 
numerical  modeling. 


PASC  numbers:  43.30. Pc,  43.30.  Hw,  43.30.Qd 


I.  INTRODUCTION 


Acoustic  ice  reflection  and  scattering  have  been  subjects  of  many  experimental  and  theoretical 
studies  since  the  earlier  1960s,  because  the  interaction  of  sound  waves  with  the  rough  sea  ice 
cover  is  the  main  factor  that  affects  attenuation  of  acoustic  signals  propagated  in  the  Arctic 
Ocean.  The  signal  transmission  loss  increases  rapidly  with  frequency  due  to  ice  scattering. 
Therefore  long-range  acoustic  transmissions  in  the  Arctic  Ocean  are  possible  only  at  low 
frequencies  of  tens  of  Hz,  and  consequently  low-frequency  acoustic  propagation  is  of  particular 
research  interest.  The  total  energy  loss  of  transmitted  signals  was  studied  in  several  acoustic 
propagation  experiments  in  the  Arctic  .  Attenuation  of  individual  normal  modes  at  low 
frequencies  is  much  less  investigated  experimentally.  Knowing  the  modal  propagation  loss  is 
important,  because  this  makes  it  possible  to  predict  the  range  and  depth  dependence  of  the  sound 
field.  Moreover,  theoretical  models  of  acoustic  ice  reflection  and  scattering  at  low  frequencies 
can  be  examined  in  more  detail  by  comparing  the  results  of  numerical  prediction  for  modal 
attenuation  due  to  ice  scattering  with  the  modal  propagation  loss  measured  experimentally. 

Propagation  losses  of  individual  modes  in  the  signals  transmitted  over  the  entire  Arctic 
Basin  were  measured  for  the  first  time  in  the  Transarctic  Acoustic  Propagation  (TAP) 
experiment 4-5  in  1994.  The  TAP  experiment  results  provided  us  with  important  data  that  were 
used  for  verifying  the  most  recent  models  of  icc  scattering,  propagation  loss,  and  reverberation  . 
However,  the  TAP  measurements  were  conducted  for  a  relatively  short  time  period  of  a  week  in 
April,  which  did  not  allow  us  to  study  temporal  variations  of  the  modal  propagation  loss  and 
modal  attenuation  with  changes  in  the  ice  cover  characteristics  and,  particularly,  seasonal 
variations  of  the  mean  ice  thickness  and  roughness  along  the  path.  If  this  correlation  is 


experimentally  examined  and  determined,  it  would  prove  the  capability  of  acoustics  to  remotely 
monitor  the  sea  ice  thickness  year  round  over  large  regions  of  the  Arctic  Ocean. 

Long-term  acoustic  transmissions  on  a  stationary  trans-Arctic  path  were  conducted  for 
the  first  time  in  the  Arctic  Climate  Observations  using  Underwater  Sound  (ACOUS)  experiment 
in  1998-1999.  Acoustic  observations  on  the  path  of  1250  km  long  from  the  Franz  Victoria  Strait 
to  the  Lincoln  Sea  lasted  for  14  months.  Two  of  the  ACOUS  signals  were  also  received  on  a 
vertical  array  at  ice  camp  APLIS  drifting  in  the  northern  part  of  the  Chukchi  Sea  2750  km  away 
from  the  source.  Signal  reception  on  vertical  arrays  allowed  us  to  separate  individual  modes  of 
the  signals  and  analyze  different  parameters  of  the  modes,  including  the  modal  propagation 
losses.  The  main  goal  of  the  ACOUS  experiment  was  to  study  long-term  changes  in  the  ocean 
temperature  and  heat  content  along  the  trans-Arctic  path  by  means  of  acoustic  thermometry.  In 
addition,  the  obtained  data  gave  us  an  opportunity  to  investigate  temporal  variations  of  modal 
attenuation  due  to  various  factors  including  seasonal  change  in  the  sea  ice  characteristics, 
variations  of  the  sound  speed,  and  the  mode  coupling  effects. 

It  is  widely  believed  that  scattering  of  acoustic  energy  from  the  rough  ice  cover  is  the 
only  significant  phenomenon  that  influences  acoustic  transmission  losses  and  modal  attenuation 
in  particular  in  the  deep  Arctic  Basin.  When  modeling  acoustic  transmission  losses  due  to  ice 
scattering,  researchers  commonly  use  certain  sound  speed  profiles  typical  for  the  Aretic  Basin. 
Although  the  Aretic  sound  speed  profiles  are  generally  upward  refracting  and  look  nearly  alike  in 
every  region  of  the  deep  Arctic  Basin,  actually  the  important  characteristics  of  the  sound  speed 
profile,  such  as  the  thickness  and  sharpness  of  the  upper  under-ice  acoustic  duct  above  the 
thermocline,  vary  considerably  from  one  region  to  another  one.  Relatively  small  spatial  and 
temporal  variations  of  the  sound  speed  profile  may  lead  to  significant  change  ot  the  modal  group 


velocities  and  attenuation.  The  other  factor  that  may  affect  strongly  the  modal  propagation  loss 
and,  in  certain  cases,  the  total  acoustic  transmission  loss  is  mode  coupling  over  range-dependent 
sections  of  the  acoustic  path.  The  influence  of  these  factors  on  the  modal  propagation  loss  on  the 
ACOUS  path  is  considered  in  this  paper  along  with  the  analysis  of  the  temporal  variations  of 
modal  attenuation  due  to  seasonal  change  in  the  ice  characteristics.  The  results  of  numerical 
modeling  arc  used  to  interpret  the  experimental  observations. 

II.  ACOUS  EXPERIMENT 
A.  Experimental  scheme 

The  ACOUS  experiment  started  in  October  1998  with  the  installation  of  an  autonomous  acoustic 
source  deployed  at  the  northern  edge  of  the  Franz  Victoria  Strait  approximately  200  km  north¬ 
west  of  Franz  Josef  Land  and  an  autonomous  vertical  receiving  array  deployed  in  the  Lincoln 
Sea  1250  km  away  (Fig.l).  Both  the  source  and  the  receiving  array  were  installed  at  the  edge  of 
the  continental  shelf.  The  source  was  moored  at  a  depth  of  60  m  below  the  sea  surface  in  water 
depth  of  440  m  and  transmitted  tomographic  signals  at  the  carrier  frequency  of  20.5  Hz  with  a 
±60°  phase  modulation  by  10  periods  of  a  255-bit  maximum  length  pseudo-random  sequence  (M- 
sequence).  Every  bit  of  the  sequence  consisted  of  10  periods  of  the  carrier  frequency.  The 
transmissions  were  scheduled  to  start  at  00:00  GMT  every  four  days.  The  total  duration  of  the 
signals  was  about  20.7  minutes  and  the  signal  level  was  nominally  195  dB  re.  1  pPa  (effective)  at 
I  m.  Two  accelerometers  installed  on  the  membranes  of  the  source  measured  the  amplitude  of 
oscillations  at  the  very  beginning  of  each  transmission  and  then  the  emitting  system  controller 
adjusted  the  oscillation  amplitude  to  the  nominal  level  of  the  transmitted  signal.  The  octal  code 
of  the  modulating  M-sequenee  is  703.  if  the  adjusted  amplitude  corresponds  to  the  nominal  level 


within  195-195.5  dB.  Otherwise  the  signal  should  be  modulated  with  a  different  M-sequence,  of 
which  the  octal  code  is  set  according  to  the  level-to-code  reference  table. 

The  acoustic  signals  from  the  source  were  received  on  the  vertical  array.  This  flexible  array 
consisted  of  8  hydrophones  spaced  at  70-m  intervals  starting  at  12  m  from  the  bottom  and  nearly 
filling  the  water  column  of  545  m  deep.  The  array  was  also  equipped  with  5  self-recording 
micro-CTD's  that  recorded  pressure,  temperature,  and  salinity  at  five  different  depths  every  10 
minutes.  In  March  2001  the  Lincoln  Sea  receiving  array  was  recovered.  Examination  of  the 
signal  recordings  on  the  array  showed  that  the  source  operated  suitably  for  14  months  from 
October  11,  1998  until  December  8,  1999.  For  this  time  period  the  source  transmitted  a  total  of 
107  signals.  The  M-scqucnce  code  was  703  in  all  of  the  transmitted  signals,  which  meant  that  the 
level  of  the  transmitted  signal  remained  nearly  nominal  (195  dB)  for  all  the  period  of  acoustic 
observations. 

On  April  9  and  13,  1999,  two  recordings  of  the  ACOUS  transmissions  were  made  on  a 
vertical  array  at  ice  camp  APLIS  established  at  the  northern  edge  of  the  Chukchi  Sea  to  support 
the  SCICEX-1999  expedition.  This  flexible  array  was  suspended  from  sea  ice  and  consisted  of  8 
equidistant  hydrophones  that  spanned  the  water  column  from  150  to  675  m.  The  array  tilt  was 
tracked  via  an  ice-mounted  acoustic  navigation  system. 


B.  Environmental  conditions 

The  bathymetry  and  the  typical  sound  speed  profile  along  the  acoustic  path  to  the  Lincoln  Sea 
are  shown  in  Fig.2.  The  path  crosses  the  Eurasia  continental  slope,  the  Gakkcl  Ridge,  the  western 
edge  of  the  Lomonosov  Ridge  and  the  Canadian  continental  slope.  Both  continental  slopes  arc 
steep  so  that  eoupling  of  the  acoustic  modes  propagated  over  the  slopes  is  significant.  The 


Gakkel  Ridge  is  not  high  enough  to  influence  the  low-order  modes  of  a  20-Hz  signal  propagated 
in  the  Arctic  upward  refracting  acoustic  channel.  The  western  edge  of  the  Lomonosov  Ridge 
rises  up  to  a  depth  of  1300  m  and  influences  some  of  the  acoustic  modes  that  carry  a 
considerable  portion  of  the  energy  of  the  propagated  signals. 

According  to  the  recent  concept  of  the  Atlantic  Intermediate  Water  (A1W)  circulation  in 

the  Arctic  Ocean 7  shown  in  Fig.l,  the  acoustic  path  crosses  the  major  circulation  gyre  along  the 
continental  margins  in  the  regions  of  the  incoming  flow  north  of  the  Barents  Sea  and  the  return 
flow  north  of  the  Canadian  shelf.  The  path  also  crosses  the  secondary  branches  of  AIW 
circulation  along  the  Lomonosov  and  Gakkel  Ridges.  The  AIW  circulation  along  the  major 
cyclonic  gyre  and  the  subsidiary  gyres  is  assumed  to  occur  with  so-called  boundary  currents 

along  the  continental  shelf  and  the  sea  ridges  ®.  The  core  of  the  AIW  flow  along  the  Eurasia 
continental  slope  is  clearly  seen  in  the  temperature  section  built  along  the  acoustic  path  with  the 
use  of  the  available  oceanographic  data.  A  statistical  analysis  of  the  historical  oceanographic  data 
collected  in  the  regions  along  the  acoustic  path  has  shown  that  the  AIW  temperature  is  most 
variable  along  the  initial  section  of  the  acoustic  path  in  the  Nansen  Basin,  and,  particularly,  over 
the  Eurasia  continental  slope.  It  was  found  that  Atlantic  waters  in  this  region  in  the  1990s  were 
noticeably  warmer  than  that  in  1970-80s.  In  the  Fram  Basin  and  the  Lincoln  Sea  the  temperature 

9 

profiles  arc  more  stable  and  vary  a  small  amount  relative  to  the  climatology  profiles  . 

The  DMSP  SSM/I  data  on  ice  concentration'0  show  that  during  the  experimental 
observations  the  sea  surface  was  covered  year-round  with  ice  in  every  region  along  the  acoustic 
path.  The  mean  ice  draft;  was  expected  to  be  approximately  2  m  in  the  Nansen  Basins  in  the  so- 
called  Arctic  summer  season  (July-November)  .  The  other  regions  along  the  path  and  the 


Arctic  winter  season  (February  -  May)  are  not  well  represented  with  the  recent  submarine  up- 
looking  sonar  data. 

The  acoustic  path  to  APLIS  crosses  the  Eurasia  continental  slope,  the  Lomonosov  and 
Mendeleyev  Ridges,  and  the  Chukchi  Plateau  (Fig.3).  The  Lomonosov  and  Mendeleyev  Ridges 
do  not  rise  above  a  depth  of  1000  m  and  therefore  do  not  influence  the  lower  modes  of  the 
ACOUS  signal.  The  sea  floor  topography  in  the  region  of  the  Chukchi  Plateau  is  very  uneven. 
The  acoustic  path  crossed  some  parts  of  the  plateau  over  which  the  depth  was  less  than  600  m 
and  hence  the  sea  floor  influenced  considerably  every  mode  of  the  ACOUS  signal  except  for 
mode  1.  The  Scientific  Ice  Expeditions  (SCICEX)  oceanographic  transect  conducted  by  USS 
Hawkbill  in  April  1999  was  almost  coincident  with  the  ACOUS  path  to  APLIS  and  made  right 
after  the  ACOUS  signal  receptions  at  APLIS.  Therefore  the  sound  speed  profiles  were  well 
defined  everywhere  along  the  acoustic  path  except  for  the  initial  400-km  section  that  was  not 
spanned  by  the  oceanographic  section.  The  oceanographic  data  collected  in  the  regions  along  this 
section  in  the  1990s  and  the  SCICEX- 1999  data  were  used  to  synthesize  the  sound  speed  section 
along  the  acoustic  path  (shown  in  Fig.3)  for  numerical  modeling  of  the  ACOUS  signal 
propagation.  In  April  1999,  every  region  along  the  ACOUS  path  to  APLIS  was  almost  entirely 

covered  with  sea  ice.  Since  the  upward  looking  sonar  data  on  ice  draft  collected  in  the  SCICEX- 
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1999  expedition  were  unavailable  for  us,  we  used  the  data  of  the  earlier  submarine  cruises  to 
estimate  the  ice  statistics  parameters  for  modeling  acoustic  propagation  along  the  path. 


III.  MODEL  OF  THE  ACOUS  SIGNAL  PROPAGATION 
A.  Ice  scattering  loss 

The  coherent  reflection  coefficient  from  Arctic  sea  ice  was  computed  with  the  use  of 
Kudrvshov's  ice  scattering  model  .  of  which  a  practical  implementation  is  discussed  in  Ref.  15. 


Proper  determination  of  the  ice  statistic  characteristics  including  the  roughness  spectrum  is 
critical  for  icc  reflection  modeling  '6.  The  morphological  characteristics  of  the  Arctic  ice  canopy 
have  changed  substantially  over  the  last  2  decades  Therefore  we  analyzed  the  data  on  Arctic 
ice  draft  which  were  available  from  the  most  recent  submarine  upward-looking  sonar 
observations,  in  order  to  determine  the  most  likely  spectrum  of  the  ice  roughness  and  the  typical 
numbers  of  mean  ice  thickness,  standard  deviation,  and  roughness  correlation  length.  The 
analysis  showed  that  the  spectral  fall-off  rate  of  2  is  most  probable  for  a  1-D  ice  draft  profile  in 
the  Arctic  winter  season,  while  it  is  somewhat  smaller  and  less  definite  in  the  summer  season 

(Fig.4).  For  the  Goff-Jordan  model  of  the  roughness  spectrum  ,  the  frequency  dependence/  of 
a  1-D  roughness  spectrum  corresponds  to  a  fractal  dimension  of  2.5  and  the  normalized  2-D 
power  spectrum 

G(u)=(2rt)_2/02[l  +  (u/0)2]  V2 ,  (0 

where  /„  is  the  roughness  correlation  length.  The  roughness  power  spectrum  (1)  corresponds  to 
azimuthal  isotropic  roughness  with  the  correlation  function  /?(/)  =  cxp(-|r|//0).  Figure  5  shows 

an  example  of  the  measured  power  spectrum  of  icc  draft  roughness  and  its  approximation  by  the 
model  power  spectrum  that  corresponds  to  the  exponential  correlation  function. 

The  distribution  of  the  ice  roughness  correlation  length,  derived  from  the  icc  draft  data,  is 
demonstrated  in  Fig.6.  In  the  winter  season,  the  ice  roughness  correlation  length  of  35-40  m  is 
apparently  most  likely.  In  Arctic  summer,  the  correlation  length  of  about  40  m  is  also  frequent, 
but  the  variance  is  much  greater  than  that  in  winter.  A  large  variance  of  the  icc  roughness 
correlation  length  is  the  consequence  of  ice-free  zones  spreading  in  summer. 


The  ice  statistic  parameters,  such  as  mean  ice  draft  (mean  thickness)  and  roughness 
standard  deviation,  vary  considerably  over  the  Arctic  Ocean  regions  and  the  seasons.  According 
to  the  most  recent  observations,  in  winter  the  mean  ice  thickness  is  about  2-3  m  in  the  Canada 
Basin  and  the  Beaufort  Sea  and  3  -  4  m  in  the  Fram  and  Nansen  Basins.  In  summer,  the  mean  ice 
thickness  decreases  to  1.5-2  m  in  the  Canada  Basin  and  2-3  m  in  the  Fram  and  Nansen  Basins.  In 
contrast  to  the  ice  thickness,  the  relationship  between  the  roughness  standard  deviation  and  the 
mean  thickness  of  ice  is  much  less  variable.  In  winter,  the  standard  deviation  is  most  likely  0.6  - 
0.65  of  the  mean  thickness,  which  is  shown  in  Fig. 7(a).  In  summer  the  values  from  0.5  to  0.7  are 
almost  equiprobable  (Fig.7(b)). 

Roughness  of  the  top  surface  of  the  Arctic  ice  cover  is  also  taken  into  account  in 
Kudryashov’s  ice  scattering  model.  Mutual  statistical  characteristics  of  the  top  and  bottom  ice 

,  1  16’18  u  u 

roughness  arc  poorly  investigated.  It  is  assumed  in  some  recent  ice  scattering  models  that  the 
top  and  bottom  scattering  surfaces  of  Arctic  ice  arc  uncorrclatcd  with  each  other.  However, 
according  to  the  Archimedean  principle,  the  roughness  of  both  surfaces  of  a  floating  icc  sheet 
must  be  partly  correlated  and  have  similar  spectra  especially  at  smaller  wavenumbers  that 
correspond  to  the  spatial  periods  larger  than  the  icc  thickness.  Local  observations  at  drifting  icc 
stations  North  Pole  have  shown  that  the  correlation  between  the  top  and  bottom  icc  roughness  is 
about  0.7,  and  the  standard  deviation  of  the  top  surface  is  about  four  times  less  than  that  of  the 
bottom  surface.  Although  these  observations  were  quite  local  and  sparse,  and  hence,  could  not 
create  a  representative  statistical  assembly,  we  used  these  estimates  in  the  icc  scattering  model  in 
the  absence  of  more  reliable  data. 

Finurc  8(a)  demonstrates  the  magnitude  of  the  ice  reflection  coefficient  versus  grazing 
angle  calculated  at  six  frequencies,  spaced  at  4  Hz  from  16  Hz  to  36  Hz,  for  the  ice  parameters 


expected  in  winter  in  the  Central  Arctic  Basin.  Small  spikes  at  a  grazing  angle  of  60°  coincide 
with  the  critical  angle  for  the  longitudinal  waves  in  ice.  The  angular  dependence  of  the  coherent 
ice  reflection  coefficient  shown  in  Fig.8(b)  was  calculated  for  the  typical  ice  characteristics 
observed  in  the  Nansen  Basin  in  August-October  in  the  1990s.  The  modal  attenuation 
coefficients  along  the  ACOUS  paths  were  calculated  for  the  typical  sound  speed  profiles  using 
KRAKENC  code  '9  with  the  modeled  ice  reflection  coefficients  as  the  top  boundary  conditions. 
Figures  9(a)  and  9(b)  show  the  variation  of  the  attenuation  coefficients  of  modes  1-5  at  20.5  Hz 
along  the  path  to  the  Lincoln  Sea  calculated  for  the  ice  parameters  given  in  the  captions  of  Figs 
8(a)  and  8(b)  respectively.  The  total  propagation  loss  of  modes  on  the  path  are  given  in  these 
figures.  It  is  evident  from  the  comparison  of  modeling  results  in  Fig. 9,  that  seasonal  change  of 
the  ice  cover  may  lead  to  noticeable  change  in  the  modal  propagation  loss  on  a  1-Mm  transarctic 
path,  especially  for  mode  1.  Note  that  the  attenuation  coefficient  of  individual  modes  also  varies 
with  the  horizontal  variation  of  the  sound  speed  profile  along  the  path,  which  results  in  variation 
of  the  modal  cycle  lengh.  Substantial  change  of  the  modal  attenuation  due  to  sound  speed 
variations  will  be  discussed  in  more  details  in  the  next  section. 

It  is  evident  that  the  ice  statistic  parameters,  such  as  mean  thickness  and  standard  deviation, 
were  not  uniform  along  the  ACOUS  path  to  the  Lincoln  Sea.  It  is  known  from  the  earlier 
observations  of  ice  draft  in  the  1970- 1980s  that  the  mean  thickness  and  standard  deviation  of  sea 
ice  in  the  Lincoln  Sea  and  the  adjacent  regions  were  larger  than  those  in  the  Nansen  and  Fram 
Basins.  However,  wide-ranging  observations  of  ice  draft  north  of  Greenland  and  the  Canadian 
Archipelago  were  not  conducted  in  the  1990s.  Therefore,  when  modeling  the  ACOUS  signal 
propaitation.  we  assumed  that  the  ice  cover  along  the  final  300-400  km  section  of  the  acoustic 
path  might  be  somewhat  thicker  than  that  in  the  Nansen  Basin. 


B.  Effects  of  mode  coupling  and  variation  of  the  sound  speed  profile 

The  propagation  of  an  acoustic  signal  at  20  Hz  from  the  Franz  Victoria  Strait  to  the  Lincoln  was 
modeled  earlier,  before  the  experiment,  using  the  climatology  sound  speed  profiles  and  the 
ET0P05  bathymetry  database.  The  coupled-mode  algorithm  used  for  numerical  modeling  and 
the  results  of  modeling  are  thoroughly  considered  in  Ref.20,  where  it  is  shown  that  the  mode 
coupling  influences  considerably  the  modal  propagation  loss  along  the  path.  It  is  also  shown  in 
Rcf.20  that  temporal  variations  of  the  sound  speed  over  any,  rapidly  or  slowly  varying  section  of 
the  path  may  lead  to  large  variations  of  the  modal  propagation  loss  due  to  changes  in  the  mode 
coupling  effects.  The  newer  modeling  results  presented  below  were  obtained  for  the  actual 
location  of  the  ACOUS  source  and  the  Lincoln  Sea  array  with  the  use  of  the  recent  International 
Bathymetry  Chart  of  the  Arctic  Ocean  (1BCAO)  database  that  is  much  more  accurate  than 
ETOP05.  Moreover  the  real  sound  speed  profiles  used  in  the  new  model  are  much  more  realistic 
than  the  climatology  data. 

Acoustic  interaction  with  the  sea  floor  plays  a  significant  role  in  attenuation  of  the  signals 

propagated  over  shallow-water  sections  of  the  acoustic  path.  For  the  ACOUS  paths,  such 

sections  are  the  northern  edge  of  the  Franz  Victoria  Strait  with  the  Eurasia  continental  slope,  the 

northern  edge  of  the  Lincoln  Sea  with  the  Canada  continental  slope,  and  the  Chukchi  Plateau. 

Over  these  regions  almost  every  mode  of  the  ACOUS  signal  underwent  significant  interaction 

with  the  seabed.  In  the  deep  basins,  interaction  of  the  ACOUS  signals  with  the  seabed  was 
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negligible.  Acoustic  properties  of  the  seabed  in  the  Lincoln  Sea  and  Franz  Victoria  Strait  are 
known  to  some  extent,  while  for  the  region  of  Chukchi  Plateau  we  have  no  geoacoustic  data 
sufficiently  reliable  for  adequate  acoustic  modeling.  Therefore  the  geoacoustic  characteristics  of 


the  Lincoln  Sea  bottom  were  extended  in  the  model  to  the  shallow-water  regions  with  unknown 
acoustic  properties  of  the  seabed. 

Figure  10  demonstrates  the  range  dependence  of  the  amplitudes  of  modes  1-5  of  the 
ACOUS  signal  at  the  central  frequency  along  the  acoustic  path  to  the  Lincoln  Sea,  disregarding 
the  factor  of  cylindrical  spreading.  This  pattern  was  calculated  for  a  particular  sound  speed 
section  that  was  combined  from  1)  the  oceanographic  data  of  a  100-km  long  oceanographic 
section  across  the  Eurasia  continental  slope  north  of  the  Franz  Victory  Strait  in  1993,  2)  several 
temperature/salinity  (T/S)  profiles  of  the  1990s  at  distances  from  150  to  550  km,  3)  the 
climatology  T/S  profiles  from  600  to  1150  km,  and  4)  the  1CESHELF  oceanographic  data  of 

1990s  in  the  Lincoln  Sea  3.  The  figure  shows  that  every  mode  of  the  ACOUS  signal  experiences 
strong  coupling  with  the  other  modes  over  the  initial,  shallow  water  section  of  the  continental 
slope.  The  mode  coupling  reduces  rapidly  with  the  increase  of  sea  depth.  Relatively  small 
variations  of  the  temperature  profiles  near  the  edge  of  the  continental  shelf  in  the  environmental 
model  (such  as  alteration  for  the  T/S  profiles  measured  in  1998  for  the  initial  20-km  section)  lead 
to  considerable  changes  of  the  modal  amplitudes  at  the  beginning  of  the  deep-water  section, 
which  is  clearly  seen  when  comparing  the  modal  amplitudes  along  the  initial  100-km  section 
shown  in  Fig.  10  with  those  in  Fig.l  1  calculated  for  the  sound  speed  profiles  of  1998.  Note  that 
the  initial  modal  amplitudes  (modal  excitation  coefficients)  arc  not  so  variable  under  possible 
changes  in  the  sound  speed  profile  at  the  source  location.  In  deeper  water,  at  distances  from  30  to 
70-80  km  modes  1-4  remain  coupled  and  hence  varying  in  amplitude  due  to  strong  horizontal 
gradients  of  the  sound  speed  in  the  intermediate  water  layer.  The  eoupling  of  modes  1-8  that 
contribute  most  to  the  energy  of  propagated  signals  is  negligible  at  distances  from  100  to  1000 


km.  Along  the  final  section  of  30-50  km  long,  almost  every  mode,  except  mode  1,  undergoes 
considerable  amplitude  variations  due  to  mode  coupling  over  the  Canada  continental  slope. 

Mesoscale  variations  of  the  sound  speed  profile  along  the  path,  which  are  large  enough  to 
influence  the  mode  coupling  effects  over  the  range-dependent  sections,  cause  short-term 
variations  of  the  modal  propagation  loss  at  the  end  of  the  path.  In  the  presence  of  intense  mode 
coupling,  the  amplitude  variations  of  modes  are  not  correlated  with  changes  in  the  sound  speed 
and  have  a  quasi-random  character.  These  short-term,  random  fluctuations  may  be  superimposed 
on  long-term  changes  in  the  modal  propagation  loss  due  to  change  in  the  ice  statistics  parameters 
along  the  path  and  hence  the  ice  scattering  strength.  Theoretical  modeling  of  coupled  mode 
scattering  in  time  varying  environmental  conditions  is  complicated,  because  it  requires  an 
adequate  stochastic  formulation  of  the  mesoscale  variations  in  the  sound  speed  field  for  solving  a 
complex  system  of  stochastic  differential  equations  for  the  coupled-mode  propagation.  For  the 
experimental  data,  it  may  be  possible  to  decompose  the  short-term  variations  and  long-term 
changes  by  means  of  spectrum  analysis. 

As  shown  in  Ref.20,  the  greater  the  distance  from  the  source  to  the  mode  coupling  zone, 
the  wider  the  modal  pulse  will  be  on  the  receiving  array.  If  the  zone  of  intense  mode  coupling  is 
far  from  the  source  and  the  group  velocities  of  coupled  modes  arc  substantially  different,  like 
those  of  modes  1  and  2  at  20  I  Iz,  the  modal  arrivals  arc  split  into  multi-peak  pulses.  In  the  case 
of  mode  coupling,  the  pulse  of  a  separate  mode  perfectly  filtered  on  the  array  actually  carries  the 
energy  of  several  modes  that  exchange  their  energy  when  propagating  over  the  range-dependent 
sections.  These  modes  may  have  different  attenuation  due  to  ice  scattering  over  the  major  portion 
of  the  path.  However,  the  resulting  propagation  loss  of  a  particular  mode  at  the  receiving  site  will 
not  be  strictly  proportional  to  the  attenuation  coefficient  of  this  mode  along  the  path.  This  makes 


it  more  difficult  to  interpret  the  results  of  the  propagation  loss  measurements  for  individual 
modes. 

In  the  oceanographic  conditions  typical  for  the  deep-water  Arctic  Basin,  mode  1  at  20  Hz  is 
trapped  in  the  low-speed  upper  Arctic  acoustic  duct  that  is  formed  by  the  sea  surface  and  the 
thermocline  between  the  upper  layer  of  cold  Arctic  water  and  the  AIW  layer.  As  a  result,  this 
mode  propagates  considerably  slower  and  attenuates  faster  than  the  higher,  mid-depth  refracted 
modes.  However,  under  certain  environmental  conditions  feasible  in  the  regions  along  the 
Eurasia  continental  slope  north  of  Spitsbergen  and  Franz  Josef  Land,  when  the  maximum  AIW 
temperature  exceeds  2.5°  C,  and  the  AIW  layer  expands  above  a  100-m  depth,  the  upper  duct 
cannot  trap  mode  1  at  20  Hz,  which  makes  this  mode  propagating  faster  and  attenuating  slower 
than  mode  2.  From  the  modal  travel  times  and  the  width  of  arrival  pulses  measured  in  the 
ACOUS  signals  on  the  Lincoln  Sea  array,  it  was  concluded  that  such  atypically  warm  conditions 
arose  on  the  ACOUS  path  over  a  wide  zone  250  -  300  km  north  of  the  source  in  the  last  quarter 

of  1 999 24.  Figure  12  shows  the  attenuation  coefficients  of  modes  1-5  versus  frequency  calculated 
for  the  typical  (a)  and  “warm”  (b)  temperature  and  sound  speed  profiles  observed  in  the  Nansen 
Basin  north  of  the  ACOUS  source  (Ref.  24,  Fig.6).  The  effective  frequency  band  of  the  ACOUS 
signal  is  approximately  20.5±1  Hz.  Hence  relatively  small  and  likely  variations  of  the  AIW 
characteristics  lead  to  change  of  the  critical  frequency  of  modes  1  in  the  upper  acoustic  duct, 
which  greatly  affects  the  attenuation  coefficient  of  this  mode  of  the  ACOUS  signal.  In  contrast  to 
mode  1,  the  attenuation  of  mode  2  increases  in  this  transient  frequency  band.  In  the  acoustic 
propagation  model  with  the  typical  sound  speed  profiles  along  the  path  (Fig. 2),  such  inversion  in 
the  attenuation  rate  of  modes  1  and  2  is  only  seen  within  a  short  interval  at  the  very  beginning  of 
the  path  (Figs  10  and  11)  and  does  not  greatly  influence  the  total  propagation  loss  of  these 


modes.  However,  if  the  warm  profiles  extend  over  300  km,  the  total  propagation  loss  of  mode  1 
decreases  by  3-4  dB  while  the  propagation  loss  of  mode  2  increases  by  2-3  dB. 

Variations  of  the  relative  amplitudes  of  modes  1-5  modeled  along  the  acoustic  path  to  ice 
camp  APLIS  are  shown  in  Fig.  13.  Mode  coupling  is  significant  for  every  mode  along  the  initial 
section  over  the  continental  slope.  The  Lomonosov  Ridge  (at  950-1000  km)  does  not  influence 
mode  1,  affects  slightly  modes  2  and  3,  and  causes  strong  coupling  of  modes  4,  5  and  higher.  For 
the  accepted  model  of  bathymetry  along  the  path,  the  Chukchi  Rise  should  not  influence  mode  1. 
The  other  modes  are  subject  to  significant  coupling  over  the  Chukchi  Rise. 


IV.  EXPERIMENTAL  RESULTS 
A.  Processing  of  the  ACOUS  signals 

After  signal  demodulation  and  pulse  compression,  the  modal  propagation  loss  can  be  estimated 
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in  the  adiabatic  approximation  by  : 
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where  pm(0  is  the  pulse-compressed  waveform  of  individual  modes  filtered  on  the  array, 
max{|/?„,| }  is  the  peak  magnitude  of  the  modal  pulse,  B  =  -J 2cxp(l95/20)  is  the  peak  signal 

level  of  the  ACOUS  source,  0  =  7t/3  is  the  phase  deviation  of  the  modulation,  and 
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arc  the  modal  excitation  coefficients  that  arc  determined  from  the  modal  eigenfunctions 
Zm(0.r0,co)  and  wavenumbers  kni  (0.  co)  at  the  source  location. 


Expression  (2)  is  exact  for  an  ideal  triangular  pulse  resulting  after  pulse  compression  of  the 
phase-coded  signal  in  an  unlimited  frequency  band.  Actually,  the  pulse-compressed  modal 
arrivals  filtered  from  the  experimental  signals  cannot  be  absolutely  triangular,  because  the 
frequency  bandwidth  of  the  source-to-receiver  transfer  functions  is  limited.  This  reduces  by  a 
considerable  amount  (depending  on  the  frequency  band  of  signal  processing)  the  peak  magnitude 
of  the  arrival  pulse.  Moreover,  the  frequency  dependence  of  the  modal  wavenumbers  is  not 
absolutely  linear,  which  may  cause  a  widening  of  the  modal  pulse  and  a  reduction  of  its 
maximum.  The  energy  of  the  modal  arrival  pulse  is  invariant  under  changes  in  the  pulse  shape 
and  hence  can  be  used  as  a  more  accurate  measure  of  the  modal  propagation  loss  than  the 
maximum  amplitude  of  the  pulse.  To  allow  for  pulse  distortion  in  equation  (2),  we  use  the 
effective  modal  amplitude  p,„  instead  of  max{|pm| }.  This  effective  amplitude  is  defined  from 

the  pulse  energy 
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where  the  modal  pulse  />„,  is  integrated  over  its  whole  duration  and  T  is  the  length  of  one  chip, 
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i.e.  one  digit  of  the  255-digit  M-sequcncc,  or  10  carrier  periods  in  our  case  .  Strong  mode 
coupling  also  leads  to  a  distortion  of  the  modal  pulse.  In  that  case,  the  estimate  of  the  effective 
modal  amplitude  pm  of  mode  m  includes  the  energy  contribution  from  the  other  modes 


interacted  with  this  mode  on  the  propagation  path. 


The  parameters  of  the  Lincoln  Sea  array  were  nearly  optimum  for  accurate  spatial  filtering 
of  the  8  lower  modes  that  contributed  most  to  the  sound  field  on  the  array.  The  method  ot  mode 
filtering  on  the  Lincoln  Sea  array  is  thoroughly  discussed  in  Ref.  24,  where  the  filtering 
algorithm  is  examined  with  the  use  of  the  signals  numerically  modeled  on  the  hydrophones  of 
the  array.  In  contrast  to  the  Lincoln  Sea  array,  the  vertical,  8-element  array  deployed  at  APLIS 
spanned  only  a  relatively  small  section  of  the  deep-water  Arctic  acoustic  channel.  The  vertical 
locations  of  the  array  hydrophones  were  such  that  the  matrix  of  modal  eigenfunctions  was  nearly 
singular  with  respect  to  four  of  the  eight  eigenvalues27.  Therefore  the  eigenvector  method-8  was 
used  for  filtering  the  modes  on  the  APLIS  array.  Numerical  modeling  showed  that  the  results  of 
mode  filtering  should  be  sufficiently  accurate  if  the  contribution  of  modes  higher  than  4  to  the 
sound  field  is  relatively  small.  Processing  of  the  ACOUS  signals  has  confirmed  this  assumption, 
which  will  be  shown  below. 

B.  Signals  on  the  Lincoln  Sea  array. 

All  of  the  107  ACOUS  signals  received  on  the  Lincoln  Sea  array  were  processed  by  filtering 
individual  modes.  The  modal  amplitudes  (energy),  travel  times,  and  phase  were  measured.  The 
temporal  coherence  of  the  received  signals  was  high  so  that  the  total  gain  of  signal  processing 
approached  its  maximum  of  34  dB.  The  resulting  signal-to-noisc  ratio  in  the  post-processed 
signals  on  individual  hydrophones  exceeded  40  dB  for  almost  all  of  the  ACOUS  signal 
recordings.  The  tilt  of  the  array  under  variable  currents  was  examined  using  the  data  from  the 
pressure  sensor  of  the  miero-CTD  deployed  at  the  top  of  the  array  and  a  hydrodynamic  model  of 
the  array  mooring  system.  From  the  results  of  numerical  modeling,  it  was  concluded  that  the 
array  tilt  was  negligible  for  mode  filtering  if  the  vertical  displacement  of  the  top  sensor  relative 


to  its  uppermost  position  did  not  exceed  a  threshold  of  0.5  m,  which  corresponds  to  horizontal 
displacement  of  approximately  20  m  for  the  top  of  the  Lincoln  Sea  array.  About  sixty  of  107 
ACOUS  receptions  meet  this  condition. 

We  first  examined  the  variation  of  the  total  energy  of  the  8  modes  filtered  on  the  array, 
which  is  almost  equivalent  to  the  depth-integrated  signal  energy.  All  of  the  107  received  signals 
were  used  for  this  analysis,  because  cross-talk  of  the  modes  on  the  tilted  array  influenced  the 
total  mode  energy  insignificantly.  Figure  14(a)  demonstrates  the  variation  in  the  level  of  the  total 
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mode  energy  p2  =  for  14  months  of  observations.  In  the  temporal  variation  of  the  modal 
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energy  one  can  notice  long-term  change  in  the  presence  of  superimposed  short-term  fluctuations. 
To  separate  the  long-term  change  from  short-term  fluctuations,  we  applied  the  wavelet 
dc-noising  procedure,  decomposing  the  signal  by  the  “coiflcts”  family  wavelets.  The  dashed  line 
in  Fig.  14(a)  shows  the  results  of  such  filtering.  The  filtered  long-term  change  has  an  evident 
seasonal  constituent  with  the  minimum  in  March-April  and  the  maximum  in  Scptcmbcr-Octobcr 
and  correlates  very  well  with  the  annual  cycle  of  changes  in  the  mean  sea  ice  thickness 
(Fig.  14(b))  derived  from  the  historical  data  on  ice  draft  and  the  Arctic  ice-ocean  model  .  The 
residual  short-term  variations  behave  as  a  random,  dclta-corrclatcd  process  with  a  zero  mean  and 
a  standard  deviation  of  approximately  1  dB.  These  short-term  fluctuations  arc  relatively  small, 
because  the  total  mode  energy  is  much  less  sensitive  to  changes  in  mode  coupling  due  to 
mesoscale  variations  of  the  sound  speed  profiles  than  the  energy  of  individual  modes. 

Figure  15  shows  the  gray-scale  intensity  plot  of  the  modal  arrival  pulses  (modes  1-3)  as  a 
function  of  the  travel  time  and  day  of  measurement.  The  amplitude,  energy,  and  shape  ol  the 
modal  pulses  underwent  considerable  variations  during  the  period  of  observations.  The 
short-term  variations  (from  one  transmission  to  the  next  one)  were,  most  likely,  the  result  ot 


sound  speed  variations  at  the  beginning  of  the  acoustic  path  over  the  shallow-water  section  of  the 
continental  slope,  which  influenced  strongly  the  energy  exchange  of  coupled  modes.  Mode  3  was 
less  variable  than  the  other  modes.  The  amplitude  of  mode  1  experienced  strong  short-term 
fluctuations  superimposed  on  noticeable  long-term  variations  due  to  seasonal  changes  of  the  ice 
characteristics  along  the  path.  The  main  arrival  pulse  of  mode  1  disappeared  completely  under 
noise  in  some  of  the  ACOUS  signal  receptions,  which  happened  mainly  during  the  Arctic  winter 
period  from  February  to  June  (days  80  -  260  in  Fig.  15). 

Figure  15  also  shows  that  the  travel  time  of  mode  1  was  rapidly  decreasing  for  the  last  four 
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months  of  observations  from  8.0  s  to  6.5  s  relative  to  the  arrival  of  mode  3.  It  was  concluded 
that  such  rapid  change  was  most  likely  a  consequence  of  warmer  and  shallower  Atlantic  water 
masses  crossing  the  acoustic  path  in  the  Nansen  Basin  within  the  initial  300-km  section  of  the 
acoustic  path  north  of  the  Eurasia  continental  slope,  which  was  happening  in  August-Dcccmbcr 
1999.  As  discussed  in  Section  III  (b),  the  group  velocity  of  mode  1  increases  considerably  under 
the  atypical ly  “warm”  conditions,  while  the  attenuation  of  this  mode  decreases.  This  is  an 
additional  factor  that  might  have  caused  a  3-4-dB  decrease  of  mode  1  attenuation  over  this 
section,  and  consequently  a  noticeable  increase  of  the  modal  amplitude  on  the  Lincoln  Sea  array 
in  August-Dcccmbcr  1999. 

Variations  of  the  effective  amplitude  pm  of  modes  1-3  arc  shown  in  Fig.  16.  Here  the 
samplings  do  not  include  the  events  of  excessively  large  tilt  of  the  array  at  the  time  of  signal 
reception.  Since  the  main  arrival  pulse  of  mode  1  was  greatly  fluctuating,  spread,  and  not 
definitely  detected  in  some  of  the  signal  receptions,  the  energy  of  this  mode  was  formally 
determined  as  energy  of  the  mode-filtered  signal  within  a  time  window  from  5.5  s  to  9  s  after  the 
arrival  of  mode  3.  The  dashed  line  gives  the  typical  noise  level  for  the  filtration  channel  of 


mode  1 .  The  procedure  of  wavelet  de-noising  cannot  be  correctly  applied  to  these  data  sets, 
because  the  samples  are  not  uniformly  spaced  in  time.  Moreover,  a  significant  number  of  mode  1 
samples  in  Fig.  16  do  not  represent  the  actual  level  of  the  modal  energy,  because  they  lie 
definitely  below  the  noise  level.  Despite  very  large  short-term  fluctuations,  it  is  evident  from 
Fig.  16  that  the  propagation  losses  of  mode  1  underwent  substantial  seasonal  variations:  the  pulse 
of  mode  1  can  be  surely  distinguished  in  approximately  50%  of  signal  recordings  made  in 
October-Dccember  1998;  in  February -June  1999,  mode  1  almost  disappeared  and  could  be 
detected  in  very  few  recordings;  in  August  1999,  the  amplitude  of  this  mode  increased  rapidly  so 
that  the  main  arrival  pulse  could  be  clearly  seen  in  almost  all  of  the  signal  recordings  until  the 
end  of  observations  in  December,  which  was  most  likely  the  result  of  superposition  of  the 
seasonal  ice  variations  along  the  whole  acoustic  path  and  the  rapid  change  in  the  oceanographic 
conditions  over  the  initial  250-300  km  section  occurred  in  August  1999. 

It  is  impossible  to  definitely  observe  any  long-term  constituent  in  the  temporal  variations 
of  the  propagation  loss  of  modes  2  and  3,  because  the  expected  seasonal  pcak-to-pcak  variation 
of  2-3  dB  for  those  modes  should  be  hidden  in  the  presence  of  the  short-term  fluctuations,  of 
which  the  standard  deviation  is  approximately  3  dB.  The  arrival  pulse  energy  of  modes  4-6 
(Fig.  17)  also  demonstrates  large  short-term  variations.  However,  the  seasonal  cycle  can  be 
distinguished  in  the  energy  variation  of  these  modes,  especially  modes  4  and  6.  Note  that  the 
energy  level  of  mode  5  is  noticeably  lower  on  average  than  that  of  mode  6,  while  the  modal 
spectrum  is  expected  to  be  steadily  declining  with  the  increase  of  the  mode  number  (or  modes  4 

and  higher.  The  same  peculiarity  of  the  modal  spectrum  was  observed  in  the  signals  that  were 
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received  at  ice  camp  Narwhal  in  the  Lincoln  Sea  during  the  TAP  experiment  in  1994  .  The 
frequency  band  of  the  ACOUS  signal  is  close  to  that  of  the  TAP  signal  (19.6±1.3  Hz).  Therefore 


we  assume  that  the  minimum  at  mode  5  in  both  TAP  and  ACOUS  signal  spectra  is  a 
consequence  of  a  specific  angular  dependence  of  the  sea  floor  reflection  coefficient  at  20  Hz  in 
the  Lincoln  Sea,  which  depends  on  the  elastic  properties  of  sediments  and  the  structure  of 
sediment  layers. 

Modes  2,  3,  and  4  attenuate  along  the  ACOUS  path  with  an  approximately  similar  rate. 
Therefore  we  calculated  the  total  mode  energy  and  propagation  losses  for  the  group  of  these 
modes,  which  allowed  us  to  considerably  reduce  the  short-term  fluctuations  due  to  variable 
effects  of  mode  coupling.  Figure  18  shows  variations  of  the  propagation  loss  of  the  group  of 
modes  2-4  and  separately  mode  1  on  this  path  during  14  months  of  observations.  The  seasonal 
cycle  in  the  variation  of  modes  2-4  propagation  losses  is  not  prominent,  but  visible.  In  the  Arctic 
winter  period  the  propagation  loss  of  this  modal  group  was  typically  3-5  dB,  which  is  in 
agreement  with  the  numerical  prediction  in  Fig.9(a).  In  the  summer  period  of  1999,  the 
propagation  loss  of  the  modal  group  decreased  to  1 .5-2  dB  on  average.  This  number  is  somewhat 
smaller  than  that  numerically  predicted  (Fig.9  (b))  for  the  summer  ice  parameters  given  in  the 
caption  of  Fig.8.  The  propagation  loss  of  mode  1  was  about  19-22  dB  in  winter  and  8-12  dB  in 
the  summer  period  of  1999.  The  propagation  loss  predicted  numerically  for  summer  is  about  12- 
13  dB  for  the  typical  sound  speed  profiles  along  the  path.  However,  if  the  “warm"  sound  speed 
profile  extends  over  the  initial  300-km  section  of  the  path,  the  predicted  propagation  loss  of 
mode  1  falls  to  9-10  dB. 

The  modal  propagation  loss  observed  in  winter  corresponds  to  the  path-average  attenuation 
coefficients  of  0.015-0.018  dB/km  for  mode  1  and  0.003-0.004  dB/km  for  modes  2-4.  These 
numbers  arc  similar  to  the  modal  attenuation  coefficients  measured  in  the  TAP  experiment  on 
the  Turpan-SIMI  path,  but  much  smaller  than  those  measured  on  the  Turpan-Nanvhal  path  which 


was  almost  coincident  with  the  1000-km  section  of  the  ACOUS  path  from  the  Nansen  Basin  to 
the  Lincoln  Sea.  The  integral  propagation  loss  on  the  path  to  Narwhal  was  about  30  dB  for  mode 
1  and  10-13  dB  for  modes  2-4.  The  higher  signal  attenuation  on  this  path  was  attributed  to  the 
fact  that  the  ice  cover  north  of  Greenland  and  the  Canadian  Archipelago  was  generally  thicker 
and  rougher  than  that  in  the  Central  Arctic  Basin.  The  ACOUS  results  indicate  that  the  winter  ice 
cover  has  become  considerably  thinner  on  average  in  the  Nansen  and  Fram  Basins  and  the 
Lincoln  Sea  in  1999  relative  to  1994. 

It  is  necessary  to  note  that  the  measured  propagation  loss  of  modes  2-4,  which  are  small  on 
this  path,  may  have  considerable  relative  errors  due  to  imperfect  calibration  of  the  array 
hydrophones  and  some  uncertainty  in  the  actual  level  of  the  transmitted  signal.  Indeed,  the 
calibration  accuracy  of  the  hydrophones  of  the  Lincoln  Sea  array  is  about  ±1.5  dB.  Although  the 
amplitude  of  oscillations  of  the  source  membrane  was  automatically  adjusted  to  the  reference 
sound  pressure  with  a  resolution  of  0.5  dB,  the  actual  uncertainty  of  the  absolute  source  level  is 
expected  to  be  at  least  ±1  dB.  Therefore  the  total  instrumental  error  of  the  absolute  propagation 
loss  measurements  may  amount  to  2-5  dB,  especially  as  some  of  the  instrumental  errors  contain 
regular  components. 

C.  Signals  on  the  APLIS  array 

The  arrival  patterns  of  modes  1-4  of  the  ACOUS  signal  filtered  on  the  vertical  array  at  ice  camp 
APLIS  on  9  and  13  April  1999  arc  shown  in  Figs.  19  (a)  and  (b)  respectively.  The  amplitude  of 
mode  I  shown  in  this  figure  is  multiplied  by  a  factor  of  2.  Mode  2  is  the  dominant  mode  in  the 
ACOUS  signals  received  at  APLIS.  The  arrival  pulse  of  this  mode  in  the  first  signal  is  split, 
which  testifies  to  strong  mode  coupling.  As  follows  from  the  width  of  the  split  pulse  and  the 


results  of  numerical  modeling,  strong  coupling  of  mode  2  with  the  higher  modes  occurred  along 
the  final  section  of  the  propagation  path  over  the  Chukchi  Plateau.  In  contrast  to  the  numerical 
prediction  shown  in  Fig.  19(c),  the  arrival  pulse  of  mode  3  in  these  signals  was  very  weak  and 
deformed  due  to  bottom  interaction  and  mode  coupling  over  the  Chukchi  Plateau.  The  arrival 
pulse  of  mode  1  is  weak,  but  distinguishable  in  both  signals.  Mode  4  is  practically  vanishing  in 
these  signals. 

The  propagation  loss  of  mode  1  measured  on  the  APLIS  array  is  41  dB  for  the  first  signal 
and  45  dB  for  the  second  signal,  which  corresponds  to  the  path-average  attenuation  coefficients 
of  0.015  dB/km  and  0.0165  dB/km  respectively.  These  numbers  are  similar  to  the  propagation 
loss  and  attenuation  coefficients  of  mode  1  measured  on  the  Turpan-SIMI  path  in  the  TAP 
experiment  (43.5±2.5  dB  and  0.0165±0.001  dB/km  respectively).  Moreover,  the  path-average 
attenuation  coefficients  of  mode  1  measured  in  April  1999  on  both  ACOUS  paths  arc 
approximately  the  same.  In  contrast  to  mode  1,  the  propagation  losses  of  mode  2  in  the  two 
ACOUS  sinnals  received  at  APLIS  arc  much  larger  than  that  measured  in  the  TAP  experiment 
on  the  Turpan-SIMI  path.  For  the  first,  adiabatic  arrival  pulse  of  mode  2,  the  propagation  loss 
was  approximately  21  dB  and  24  dB  in  these  two  ACOUS  signals,  versus  13±1.5  dB  measured  in 
the  TAP  experiment.  The  path-average  attenuation  coefficient  of  this  mode  estimated  on  the  path 
to  APLIS  is  0.008-0.009  dB/km,  which  is  approximately  two  times  higher  than  that  measured  on 
the  ACOUS  path  to  the  Lincoln  Sea.  The  excessive  attenuation  of  mode  2  on  this  path  was,  most 
likely,  the  result  of  strong  interaction  of  this  mode  with  the  sea  floor  in  the  region  of  the  Chukchi 
Plateau.  The  evidence  of  strong  mode  coupling  for  mode  2  and  lack  of  mode  3  in  the  ACOUS 
signals  received  at  APLIS  support  this  assumption.  An  obvious  discrepancy  between  the 
experimental  and  modeled  arrival  patterns  of  modes  on  the  APLIS  array  is  probably  a  result  of 


incorrect  modeling  of  the  bathymetry  and  clastic  properties  of  the  sea  floor  along  the  section  of 
the  acoustic  path  over  the  Chukchi  Plateau. 

Conclusions 

The  ACOUS  experiment  has  proven  the  capability  of  low-frequency  cross-Arctic  acoustic 
transmissions  to  observe  long-term  basin-scale  changes  in  the  Arctic  ice  cover.  The  seasonal 
cycle  in  the  variations  of  the  total  energy  losses  of  the  transmitted  signal  and  the  propagation 
losses  of  individual  modes  was  observed.  According  to  the  theory  of  ice  scattering,  coherent 
losses  of  the  low  frequency  signals  propagated  under  ice  in  the  Arctic  Ocean  depend  primarily 
on  the  height  and  correlation  length  of  the  ice  roughness.  As  proven  by  numerous  upward- 
looking  sonar  observations  of  Arctic  ice  draft,  the  height  of  ice  roughness  closely  mirrors  the 
change  in  the  mean  ice  thickness  through  a  nearly  linear  relationship.  For  ice  draft  in  excess  of  1 
m,  the  standard  deviation  of  the  bottom  ice  boundary  is  about  0. 5-0.6  of  the  mean  ice  thickness. 
Seasonal  melting  and  growing  of  sea  ice  result  basically  in  variations  of  the  mean  ice  thickness 
and  standard  deviation  of  ice  roughness,  which  influences  considerably  the  acoustic  propagation 
losses  and  hence  can  be  observed  by  means  of  long-range  acoustic  transmissions,  as  was  done  in 
the  ACOUS  experiment. 

The  results  of  propagation  loss  modeling  along  the  ACOUS  path  to  the  Lincoln  Sea,  using 
the  most  recent  ice  dralt  data  and  the  model  of  seasonal  change  of  ice  thickness,  arc  basically  in 
agreement  with  the  results  of  experimental  observations  on  this  path.  This  allows  us  to  derive  the 
correlation  between  the  modal  propagation  loss  and  ice  thickness  and  hence  gives  a  way  to  build 
an  acoustic  inversion  method  for  remote  observations  of  ice  thickness.  However,  the  accuracy  ot 
such  acoustic  monitoring  of  Arctic  sea  ice  needs  to  be  examined  in  the  future  experiments  by 


means  of  ice  draft  profiling  conducted  along  the  acoustic  path  at  the  time  of  acoustic 
observations. 

Numerical  modeling  has  shown  that  the  sound  speed  profile  influences  considerably  the 
attenuation  of  individual  modes  due  to  ice  scattering.  Recent  oceanographic  observations  and  the 
ACOUS  results  give  evidence  of  large-scale  spatial  and  temporal  variations  of  the  sound  speed 
profile  in  the  Arctic  Ocean.  Therefore  modelers  and  experimental  researchers  should  take  into 
account  the  dependence  of  acoustic  transmission  losses  on  the  sound  speed  profile  when 
developing  acoustic  propagation  models  and  planning  acoustic  experiments  in  the  Arctic. 

Attenuation  of  mode  1  at  20  Hz  is  much  more  sensitive  to  changes  of  the  ice  cover  than  the 
higher  modes.  Remote  observations  of  Arctic  sea  ice  by  means  of  low  frequency  cross-Arctic 
acoustic  transmissions  would  be  much  more  accurate,  if  mode  1  is  not  subject  to  significant 
interaction  with  the  sea  floor  and  strong  coupling  with  the  higher  modes  due  to  variable 
bathymetry  and  horizontal  changes  in  the  sound  speed  profile.  Such  propagation  conditions  for 
mode  1  can  be  achieved  through  choosing  cither  better  locations  ot  the  source  and  the  receiving 
array  or  a  more  suitable  frequency  band  of  the  signal.  For  example,  at  23  Hz  mode  1  would 
experience  somewhat  weaker  interaction  with  the  sea  floor  along  the  initial  section  of  the 
ACOUS  path  and  much  less  significant  coupling  with  the  higher  modes  over  the  frontal  zone  in 
the  Nansen  Basin  than  that  observed  at  20.5  Hz  in  the  ACOUS  experiment. 
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Figure  captions 


Fig.  1.  ACOUS  experiment  paths  in  the  Arctic  Ocean.  The  salmon-colored  lines  and  arrows 
show  the  recently  proposed  scheme  of  Atlantic  water  circulation  .  The  blue-colored  lines  and 
arrows  show  the  same  but  for  the  Pacific  water  circulation. 

Fig.  2.  Bathymetry  and  the  typical  sound  speed  profile  along  the  acoustic  path  to  the  Lincoln 
Sea. 

Fig.  3.  Bathymetry  and  the  sound  speed  profile  along  the  acoustic  path  to  APLIS. 

Fig.  4.  Histograms  of  the  fall-off  rate  of  the  ice  roughness  power  spectra  calculated  for  the  ice 
draft  profiles  that  were  observed  in  the  spring  of  1991  (a,  total  length  of  profiles  3500  km)  and 
the  autumn  of  1996  (b,  total  length  of  profiles  5350  km).  The  ice  roughness  spectrum  was 
calculated  for  every  profile  section.  The  length  of  sections  varied  from  several  km  to  50  km. 

Fig.  5.  One  of  the  ice  roughness  power  spectra  measured  in  1991  (dashed  line  and  dots)  and  its 
approximation  by  a  model  Golf-Jordan  spectrum  with  a  2-D  fractal  dimension  of  2.5  (solid  line). 

Fig.  6.  Histograms  of  the  ice  roughness  correlation  length  calculated  for  the  ice  dralt  profiles  that 
were  observed  in  the  spring  of  1991  (a)  and  the  autumn  of  1996  (b). 

Fig.  7.  Same  as  Fig.6.  but  for  the  ratio  of  roughness  standard  deviation  and  mean  ice  thickness. 

Fig.  8  Coherent  reflection  coefficient  versus  grazing  angle  at  frequencies  16,  20,  24,  28,  32,  and 
36  llz  for  Arctic  sea  ice  with  the  typical  winter  parameters  (a:  mean  thickness  =  3.0  in:  standard 
deviation  of  bottom  ice  boundary  =  1.8  in:  standard  deviation  of  top  ice  boundary  =  0.45  in: 
correlation  length  40  m)  and  the  typical  summer  parameters  (b:  mean  thickness  -  2.0  m: 


standard  deviation  of  bottom  ice  boundary  -  1.2  m;  standard  deviation  of  top  ice  boundary  -  0.3 
m;  correlation  length  -  40  m).  The  following  acoustic  properties  of  ice  were  used  in  the  model: 
compressional  wave  velocity  —  3000  m/s;  shear  wave  velocity  —  1800  m/s;  compressional  wave 
attenuation  -  0.45  dB/X;  shear  wave  attenuation  -  0.9  dB/X;  relative  density  -  0.91. 

Fig.  9.  Attenuation  coefficients  of  modes  1-5  at  20.5  Hz  along  the  ACOUS  path  and  the  resulting 
propagation  loss  (PL)  modeled  for  the  winter  (a)  and  summer  (b)  ice  conditions  corresponding  to 
the  ice  reflection  coefficient  shown  in  Figs. 8  (a)  and  (b)  respectively.  This  model  prediction  was 
obtained  using  the  sound  speed  profiles  along  the  ACOUS  path  typical  in  1990s. 

Fig.  10.  Modeled  variations  of  the  amplitudes  of  modes  1-5  at  the  central  frequency  of  the 
ACOUS  signal  along  the  path  to  the  Lincoln  Sea,  disregarding  cylindrical  spreading  (section 
from  100  to  1000  km  is  not  shown).  The  sound  speed  vertical  section  along  the  path  was 
combined  from  1)  the  oceanographic  data  of  a  100-km  long  oceanographic  section  across  the 
Eurasia  continental  slope  north  of  the  Franz  Victory  Strait  in  1993,  2)  several 
tcmpcraturc/salinity  (T/S)  profiles  of  the  1990s  at  distances  from  150  to  550  km,  3)  the 
climatology  T/S  profiles  from  600  to  1150  km,  and  4)  the  ICESHELF  oceanographic  data  of 
1990s  in  the  Lincoln  Sea. 

Fig.  11.  Variations  of  the  amplitudes  of  modes  1-5  at  the  central  frequency  along  the  initial  100- 
km  section  of  the  ACOUS  path  modeled  for  the  sound  speed  profiles  observed  in  1998  at 
distances  up  to  20  km  north  of  the  acoustic  source  and  the  sound  speed  profiles  measured  in 
1993  in  the  deep-water  part  of  this  section. 


Fig.  12.  Attenuation  coefficients  of  modes  1-5  versus  frequency  for  the  ice  reflection  coefficient 
shown  in  Fig.8(b)  and  the  typical  (a)  and  “warm”  (b)  sound  speed  profiles  shown  in  Ref.  25, 
Fig.6. 

Fig.  13.  Modeled  variations  of  the  amplitudes  of  modes  1-5  at  the  central  frequency  ot  the 
ACOUS  signal  along  the  path  to  APLIS,  disregarding  cylindrical  spreading.  The  sound  speed 
vertical  section  along  the  path  was  combined  from  1)  the  historical  oceanographic  data  of  1990s 
for  the  initial  350-km  section  of  the  path  and  2)  the  SCICEX-1999  T/S  profiles  observed  along 
the  major  section  from  350  km  to  2720  km. 

Fig.  14.  (a)  -  Variations  of  the  total  energy  of  modes  1-8  filtered  on  the  Lincoln  Sea  array.  The 
dashed  line  shows  the  result  of  wavelet  filtering  of  the  long-term  component  of  variations,  (b)  - 

Modeled  seasonal  cycle  of  ice  draft  (dashed)  and  thickness  (solid)  . 

Fig.  15.  Gray-scale  intensity  plot  of  the  modal  arrival  pulses  (modes  1-3)  as  a  function  of  the 
travel  time  and  the  day  of  measurements.  The  signals  arc  synchronized  with  the  arrival  time  of 
mode  3  shown  by  a  dashed  line. 

Fig.  16.  Variations  of  the  energy  of  modes  1  (blue),  2  (green),  and  3  (red)  filtered  on  the  Lincoln 
Sea  array  during  the  ACOUS  experiment.  The  dashed  line  shows  the  typical  noise  level  in  mode 
I  filtered  channel. 

Fig.  17.  Same  as  Fig.  16,  but  for  modes  4  (blue),  5  (green),  and  6  (red). 

Fig.  18.  Propauation  loss  of  mode  1  (blue)  and  the  group  of  modes  2-4  (red)  received  on  the 
Lincoln  Sea  array  during  the  ACOUS  experiment.  The  estimates  of  mode  1  propagation  loss 


above  the  dashed  line  arc  erroneous  and  hence  are  not  shown,  because  they  were  obtained  at  too 
low  signal-to-noise  ratio. 

Fig.  19.  The  arrival  patterns  of  modes  1-4  of  the  ACOUS  signal  filtered  on  the  vertical  array  at 
ice  camp  APL1S  on  9  April  at  range  R  =  2714.84  km  from  the  source  (a)  and  13  April  at  range  R 
=  2722.44  km  (b).  The  bottom  panel  (c)  shows  the  modeled  arrival  patterns  of  modes  1-6  of  the 
ACOUS  signal  on  the  APLIS  array  at  R  =  2722.44  km.  The  amplitude  of  mode  1  is  shown 
doubled  in  all  three  plots. 
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ABSTRACT 

The  Transarctic  Acoustic  Propagation  experiment  in  1994  revealed  integral,  basin-scale  warming  of  the  Atlantic  intermediate  water  layer  relative  to  climatology, 
which  supported  the  earlier  observations  of  warming  in  this  layer  in  certain  regions  of  the  Arctic  Ocean.  Both  experimental  and  modeling  results  have  shown  that 
the  travel  time  of  individual  modes  of  a  low-frequency  acoustic  signal  is  a  precise  indicator  of  changes  in  the  integral  Atlantic  water  temperature  along  cross-Arctic 
sections.  The  first  long-term  stationary  system  of  Arctic  acoustic  thermometry  was  experimentally  tested  for  14  months  in  1998-1999  in  the  framework  of  the  Arctic 
Climate  Observations  using  Underwater  sound  (ACOUS)  program.  Remote  acoustic  observations  on  the  cross-Arctic  path  from  Franz  Josef  Land  to  the  Lincoln  Sea 
detected  substantial  warming  of  Atlantic  waters  and  a  shoaling  of  the  thermocline  that  occurred  rapidly  in  the  central  Nansen  Basin  (83  -  840  N,  20  -  300  E)  In  the  last 
quarter  of  1999.  The  long-term  acoustic  transmissions  in  1998-1999  were  also  capable  of  detecting  seasonal  variations  of  the  mean  thickness  of  sea  ice  along  the 
cross- A  retie  path.  At  present,  an  extensive  network  of  acoustic  thermometry  paths  for  multiyear  observations  in  the  Arctic  Ocean  is  projected.  The  feasibility  of 
acoustic  halinometry,  l.e.  remote  observations  of  salinity  change  in  the  upper  layer,  is  also  examined  by  numerical  modeling. 


The  arrival  pattern  of  the  signal  received  at  Ice  camp  SlMl  in  the 
Transarctic  Acoustic  Propagation  experiments  In  April  1994  (blue), 
and  the  arrival  times  and  amplitudes  of  modes  1-4  (red  arrows) 
numerically  modeled  with  the  use  of  climatology  oceanographic 
data  (Tor  1970s).  Difference  In  the  experimental  and  modeled  travel 
limes  of  mode  2  of  approximately  2  s  Indicates  warming  of  almost 
.2  0  C  In  the  section-average  temperature  of  the  Atlantic  water 
layer  along  the  acoustic  path 


The  section-average  temperature  of  the  AW  layer  derived  from 
the  SCICEX  transect  of  different  years  and  the  climatology  data 
versus  the  travel  time  of  modes  2  (stars)  and  3  (cirdea)  of  an 
acoustic  signal  at  20.6  Hi  propagated  over  the  SCICEX  section, 
and  the  best  linear  fit  tor  the  relationship  between  these 
characteristics. 


Acoustic  observations  on  the  ACOUS  acoustic  path  to  the  Lincoln  Sea  (solid  red  line)  were 
conducted  tor  14  months  In  1999-1999.  Two  of  the  107  ACOUS  signals  were  received  at  Ice  camp 
APLIS  The  acoustic  path  to  APUS  (dashed  red  Hne)  was  dose  to  the  SCiCEX  transects  conducted 
In  1996,  1999,  1999.  and  2000.  A  notional  monitoring  grid  for  the  Arctic  Ocean  la  shown  by  black 
dashed  lines.  This  network  Is  called  AT  AM.  which  stands  for  Acoustic  Thermometry  and 
Autonomous  Monitoring.  In  addition  to  acoustic  receive  arrays,  the  ATAM  moorings  are  planned  to 
equip  with  oceanographic,  biological,  chemical,  sea  Ice.  and  seismic  sensors.  The  seashore 
Interface  In  the  Lincoln  Sea  already  exists  Such  an  Interface  In  the  Beaufort  Sea  does  not  yet  exist. 


•bM  ATAM  mom  mp 


a  Variations  of  the  total  energy  of  modes  1-9  filtered  on  the  Lincoln  Sea  array.  The  red 
line  shows  the  result  of  wavelet  filtration  of  the  long-term  component  of  variations .  b 
modeled  seasonal  cycle  of  Ice  draft  (dashed)  and  thickness  (solid)  [D.A  Rothrock 
et.al.,  GRL,v.26(23),  pp.  3499-3472],  A  seasonal  component  can  be  dearly  seen  In  the 
acoustic  data.  The  numerically  modeled  change  of  the  modal  propagation  loss  due  to 
seasonal  variations  of  the  Ice  thlkness  corresponds  to  the  observed  data. 


Sound  speed  in  sea  water  at  low  temperatures  At  freeing  point, 
the  sound  speed  varies  mainly  with  the  water  salinity. 
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The  amplitude  of  modes  1-3  arrival  pulses  as  a  function  of  the 
travel  time  and  the  day  of  measurements  Considerable  long¬ 
term  varlaOons  of  the  travel  time,  amplitude,  and  pulse  width  of 
mode  1  Indicated  large-scale  changes  In  the  Atlantic  water 
temperature  and  the  thermocline  depth  over  the  Initial  300-km 
section  of  the  acoustic  path.  The  varlattons  of  mode  2  travel  time 
due  to  those  changes  are  also  noticeable,  but  not  so  large. 


CONCLUSIONS 


The  temperature  sacSon  shown  In  the  upper  panel  was  modeled  using  the  oceanographic  data 
collected  along  the  Initial  500-km  section  of  the  ACOUS  path  In  the  1990s.  The  ACOUS  signal 
parameters  observed  In  the  beginning  of  the  experiment  are  In  agreement  with  this  profile 
The  bottom  panel  demonstrates  the  temperature  section  that  was  obtained  by  Inversion  of  the 
modal  characteristics  of  the  ACOUS  signals  received  In  late  November  ■  December  1999 


The  temperature  seettons  along  the  Initial  500-km  section  of  the  ACOUS  path  In  April,  August, 
and  December  of  1999  derived  from  the  results  of  numerical  slmulatton  using  the  MOM-2 
ocean  model  by  the  Geophysical  Fluid  Dynamics  Laboratory  [Karcher  at.  el..  JGR.  v.  109.  C2J. 
Note  that  the  model  also  shows  substantial  shoaling  of  the  thermocline  at  the  end  of  1999 


1.  The  first  detection  of  current  worming  of  Atlantic  waters  on  a  basin  scale  in  the  Arctic  Ocean  was  made  by  means  of  acoustic  thermometry  in  the  TAP 
experiment  in  1994. 

2.  Long-term  acoustic  observations  conducted  in  the  ACOUS  experiment  in  1998-1999  revealed  dramatic  change  that  occurred  in  the  water  masses  in  the  Nansen 
Basin  north  of  the  Franz  Victory  Strait  in  August-December  of  1999.  According  to  the  ACOUS  results,  the  Adantic  layer  become  considerably  thicker  and  warmer 
and  the  thermocline  rose  above  a  100-m  depth  in  the  Nansen  Basin  150-300-km  north  of  the  Eurasia  continental  slope. 

3.  The  ACOUS  experiment  has  proven  the  capability  of  low-frequency  cross-Arctic  acoustic  transmissions  to  observe  long-term  basin-scale  changes  in  the  Arctic 
Ice  cover.  The  seasonal  cycle  in  the  variations  of  the  total  energy  losses  of  the  transmitted  signal  and  the  propagation  losses  of  individual  modes  was  observed. 

4.  The  results  of  numerical  modeling  show  that  change  in  the  integral  salinity  of  the  upper  under-ice  layer  can  be  remotely  observed  by  means  of  acoustic 
transmissions  at  frequencies  of  80-100  Hz.  The  accuracy  of  such  measurements  is  a  subject  of  further  studies. 

The  work  was  funded  by  the  National  Science  Foundation  (Office  of  Polar  Programs),  the  Office  of  Naval  Research  (High  Latitude  Programs),  the  U.S.  Civilian 
Research  and  Development  Foundation  (Award  No.  RG2-240T),  RFBR  Grant  02-05-64339,  and  the  Ministry  of  Industry,  Science,  and  Technology  of  the  Russian 
Federation  HIIF1277.2003.S. 


Correlation  between  the  observed  salinity  of  the  upper  layer  and 
the  results  of  acoustic  Inversion  using  the  propagation  spaed  V 


The  mean  salinity  of  the  upper  30- m  water  layer  obtained  from 
the  SOCEX-1999  oceanographic  transect  (blue),  and  the  results 
of  acoustic  Inversion  for  salinity  using  the  travel  times  of  modes 
1,  3.  and  4  at  90  Hz  over  short  sections  along  the  transect. 


